REPORT  DOCUMENTATION  PAGE 


Form  Approved  OMB  NO.  0704-0188 


The  public  reporting  burden  for  this  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions, 
searching  existing  data  sources,  gathering  and  maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments 

regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information,  including  suggesstions  for  reducing  this  burden,  to  Washington 

Headquarters  Services,  Directorate  for  Information  Operations  and  Reports,  1215  Jefferson  Davis  Highway,  Suite  1204,  Arlington  VA,  22202-4302. 
Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  any  oenalty  for  failing  to  comply  with  a  collection  of 
information  if  it  does  not  display  a  currently  valid  OMB  control  number. 

PLEASE  DO  NOT  RETURN  YOUR  FORM  TO  THE  ABOVE  ADDRESS. 


1.  REPORT  DATE  (DD-MM-YYYY) 


2.  REPORT  TYPE 
New  Reprint 


4.  TITLE  AND  SUBTITLE 

Noiseless  non-reciprocity  in  a  parametric  active  device 


3.  DATES  COVERED  (From  -  To) 


5a.  CONTRACT  NUMBER 
W91  INF-09-1-0336 


5b.  GRANT  NUMBER 


6.  AUTHORS 

M.  H.  Devoret,  John  Clarke,  Archana  Kamal 


5c.  PROGRAM  ELEMENT  NUMBER 
411359 


5d.  PROJECT  NUMBER 


5e.  TASK  NUMBER 


5f.  WORK  UNIT  NUMBER 


7.  PERFORMING  ORGANIZATION  NAMES  AND  ADDRESSES 

University  of  California  -  Berkeley 

Sponsored  Projects  Office 

The  Regents  of  the  University  of  California 

Berkeley,  CA  94704  -5940 


9.  SPONSORING/MONITORING  AGENCY  NAME(S)  AND 
ADDRESS(ES) 

U.S.  Army  Research  Office 
P.O.Box  12211 

Research  Triangle  Park,  NC  27709-2211 


8.  PERFORMING  ORGANIZATION  REPORT 
NUMBER 


10.  SPONSOR/MONITOR'S  ACRONYM(S) 
ARO 


11.  SPONSOR/MONITOR'S  REPORT 
NUMBER(S) 

56202-PH-CSO.4 


12.  DISTRIBUTION  AVAILIBILITY  STATEMENT 
Approved  for  public  release;  distribution  is  unlimited. 


13.  SUPPLEMENTARY  NOTES 

The  views,  opinions  and/or  findings  contained  in  this  report  are  those  of  the  author(s)  and  should  not  contrued  as  an  official  Department 
of  the  Army  position,  policy  or  decision,  unless  so  designated  by  other  documentation. 


14.  ABSTRACT 

Non-reciprocal  devices  such  as  circulators  and  isolators  belong  to  an  important  class  of  microwave  components 
employed  in  applications  including  the  measurement  of  mesoscopic  circuits  at  cryogenic  temperatures.  The 
measurement  protocols  usually  involve  an  amplification  chain  that  relies  on  circulators  to  separate  input  and  output 
channels  and  to  suppress  backaction  from  different  stages  on  the  sample  under  test.  In  these  devices  the  usual 
reciprocal  symmetry  of  circuits  is  broken  by  the  phenomenon  of  Faraday  rotation  based  on  magnetic  materials  and 


15.  SUBJECT  TERMS 

Circulators,  isolators,  parmetric  modulation,  quantum  computing 


16.  SECURITY  CLASSIFICATION  OF: 

17.  LIMITATION  OF 

15.  NUMBER 

19a.  NAME  OF  RESPONSIBLE  PERSON 

a.  REPORT 

b.  ABSTRACT 

c.  THIS  PAGE 

ABSTRACT 

OF  PAGES 

John  Clarke 

UU 

UU 

UU 

UU 

19b.  TELEPHONE  NUMBER 

510-642-3069 

Standard  Form  298  (Rev  8/98) 
Prescribed  by  ANSI  Std.  Z39. 1 8 


Report  Title 


Noiseless  non-reciprocity  in  a  parametric  active  device 

ABSTRACT 

Non-reciprocal  devices  such  as  circulators  and  isolators  belong  to  an  important  class  of  microwave  components 
employed  in  applications  including  the  measurement  of  mesoscopic  circuits  at  cryogenic  temperatures.  The 
measurement  protocols  usually  involve  an  amplification  chain  that  relies  on  circulators  to  separate  input  and  output 
channels  and  to  suppress  backaction  from  different  stages  on  the  sample  under  test.  In  these  devices  the  usual  reciprocal 
symmetry  of  circuits  is  broken  by  the  phenomenon  of  Faraday  rotation  based  on  magnetic  materials  and  fields. 
However,  magnets  are  averse  to  on-chip  integration,  and  magnetic  fields  are  deleterious  to  delicate  superconducting 
devices.  Here  we  present  a  new  proposal  that  combines  two  stages  of  parametric  modulation  to  emulate  the  action  of  a 
circulator.  It  is  devoid  of  magnetic  components  and  suitable  for  on-chip  integration.  As  the  design  is  free  of  any 
dissipative  elements  and  based  on  reversible  operation,  the  device  operates  noiselessly,  giving  it  an  important 
advantage  over  other  non-reciprocal  active  devices  for  quantum  information  processing  applications. 


REPORT  DOCUMENTATION  PAGE  (SF298) 
(Continuation  Sheet) 


Continuation  for  Block  13 


ARO  Report  Number  56202.4-PH-CSQ 
Noiseless  non-reciprocity  in  a  parametric  active 


Block  13:  Supplementary  Note 

©2011  .  Published  in  Nature  Physics,  Vol.  7,  (4),  Ed.  0  (2011),  (Ed. ).  DoD  Components  reserve  a  royalty-free,  nonexclusive 
and  irrevocable  right  to  reproduce,  publish,  or  otherwise  use  the  work  for  Federal  purposes,  and  to  authroize  others  to  do  so 
(DODGARS  §32.36).  The  views,  opinions  and/or  findings  contained  in  this  report  are  those  of  the  author(s)  and  should  not  be 
construed  as  an  official  Department  of  the  Army  position,  policy  or  decision,  unless  so  designated  by  other  documentation. 


Approved  for  public  release;  distribution  is  unlimited. 


nature 

physics 


LETTERS 


PUBLISHED  ONLINE:  30  JANUARY  2011  |  DOI:  10.1038/NPHYS1893 


Noiseless  non-reciprocity  in  a  parametric 
active  device 


Archana  Kamal1,  John  Clarke2  and  M.  H.  Devoret1 ** 

Non-reciprocal  devices  such  as  circulators  and  isolators  belong 
to  an  important  class  of  microwave  components  employed  in 
applications  including  the  measurement  of  mesoscopic  circuits 
at  cryogenic  temperatures1-5.  The  measurement  protocols 
usually  involve  an  amplification  chain  that  relies  on  circulators 
to  separate  input  and  output  channels  and  to  suppress 
backaction  from  different  stages  on  the  sample  under  test. 

In  these  devices  the  usual  reciprocal  symmetry  of  circuits 
is  broken  by  the  phenomenon  of  Faraday  rotation  based  on 
magnetic  materials  and  fields6.  However,  magnets  are  averse 
to  on-chip  integration,  and  magnetic  fields  are  deleterious  to 
delicate  superconducting  devices7'8.  Here  we  present  a  new 
proposal  that  combines  two  stages  of  parametric  modulation 
to  emulate  the  action  of  a  circulator.  It  is  devoid  of  magnetic 
components  and  suitable  for  on-chip  integration.  As  the 
design  is  free  of  any  dissipative  elements  and  based  on 
reversible  operation,  the  device  operates  noiselessly,  giving 
it  an  important  advantage  over  other  non-reciprocal  active 
devices  for  quantum  information  processing  applications. 

Reciprocity  is  one  of  the  fundamental  symmetries  frequently 
encountered  in  electrical  circuits.  It  is  equivalent  to  the  more 
familiar  notion  of  the  principle  of  reversibility  in  optics,  which 
states  that  any  experiment  is  symmetric  under  an  exchange  of 
source  and  image9.  Reciprocity  can,  however,  be  violated,  for 
example,  by  the  magneto-optic  Faraday  effect6 — the  rotation  of  the 
polarization  vector  of  light  resulting  from  different  propagation 
velocities  of  left-  and  right-circularly  polarized  waves  in  the 
presence  of  an  applied  magnetic  field  B  parallel  to  the  direction  of 
propagation  (Fig.  1).  The  non-reciprocal  phenomenon  of  Faraday 
rotation  should  be  contrasted  with  the  superficially  similar,  though 
reciprocal,  effect  of  optical  activity  where  the  polarization  vector  of 
light  is  rotated  on  passage  through  a  non-centrosymmetric  (chiral) 
medium.  This  change  in  the  sense  of  rotation  of  polarization  for 
counterpropagating  waves  in  a  Faraday  medium  (as  seen  by  the 
observer  receiving  the  light)  has  led  some  physicists  to  refer  to  this 
effect  as  a  form  of  time-reversal  symmetry  breaking,  a  use  of  words 
that  we  prefer  to  avoid  here10. 

The  phenomenon  of  non-reciprocity  has  propelled  numerous 
theoretical  investigations  (ref.  11  and  references  therein);  further¬ 
more,  it  offers  immediate  practical  applications.  Recent  progress 
in  solid-state  superconducting  qubits,  which  provide  some  of  the 
most  promising  architectures  for  scalable  quantum  computers12’13, 
has  generated  a  strong  incentive  to  integrate  the  components 
required  for  qubit  operations  and  readout  on-chip  for  incorpora¬ 
tion  in  future  quantum  mechanical  processors.  A  large  variety  of 
qubit  readout  protocols  involve  microwave-reflection-based  mea¬ 
surements  and  rely  on  non-reciprocal  devices  such  as  circulators 
(or  isolators)  for  separation  of  input  and  output  channels3,5,14. 
These  devices  also  play  a  strategic  role  in  measurements  based 


Figure  1  |  Faraday  rotation  and  circulator  action,  a,  Faraday  rotation  for  a 
wave  travelling  from  left  to  right  in  a  Faraday-active  medium,  followed  by  a 
reflection  back  into  the  medium  leading  to  a  reversal  of  the  direction  of 
propagation.  The  rotation  of  the  light  polarization  is  fixed  to  a  rotation-like 
property  of  the  medium  (shown  by  the  arrows),  set  by  an  external  magnetic 
field  oriented  along  the  propagation  axis.  The  sense  of  light  rotation  as  seen 
with  respect  to  the  direction  of  propagation  remains  the  same,  leading  to 
the  doubling  of  the  rotation  angle  on  reversing  the  ray  through  the  medium, 
b,  Rotation  of  the  polarization  vector  of  light  on  passage  through  an 
optically  active  medium,  on  the  other  hand,  cancels  out  on  reversing  the 
direction  of  propagation.  This  occurs  because  optical  rotation  depends  on 
the  chirality  of  the  medium  (represented  as  a  helix),  which  also  reverses 
with  the  direction  of  propagation,  c,  Representation  and  schematic  design 
of  a  conventional  four-port  circulator.  The  device  consists  of  two  90° 
hybrids  (equivalent  to  optical  beam  splitters)  separated  by  a  non-reciprocal 
phase  shifter  based  on  Faraday  rotation.  The  filled  black  arrows  indicate  an 
amplitude  split  with  no  phase  change  and  the  open  arrows  indicate  an 
amplitude  split  with  a  90°  phase  change.  The  non-reciprocal  phase  shift  is 
effective  only  for  the  propagation  direction  indicated  by  the  arrow  on  the 
phase-shifter  rectangle. 
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Figure  2  |  Description  of  an  active  reversible  (information-conserving) 
in-phase/quadrature  modulator  carrying  out  frequency  UDC.  a,  Circuit 
schematic  of  the  UDC  containing  only  dispersive  components.  The  two 
low-frequency  series  LC  resonators  (with  L/\  =  Lb  and  C a  =  Cb)  are  fed  by 
two  input  semi-infinite  transmission  lines,  A  and  B,  and  parametrically 
coupled  to  a  third  high-frequency  series  LC  resonator  leading  to  an  output 
line  C.  The  parametric  coupling  is  achieved  by  varying  the  mutual 
inductances  Mi  and  M2  between  the  left  and  right  resonators  at  the  carrier 
frequency  coc,  which,  for  optimal  frequency  conversion,  is  set  at  the  band 
centre  of  the  right  resonator.  When  operated  from  left  to  right,  the  circuit 
carries  out  the  modulation  of  low-frequency  signals  of  frequency  u>q 
travelling  on  ports  A  and  B  to  generate  sidebands  at  a>c  ±a>o  travelling  on 
the  high-frequency  line  C.  It  carries  out  the  inverse  operation  of 
demodulation  when  operated  in  reverse  from  right  to  left,  b,  Spectral 
density/response  landscape  for  different  spatial  channels  (or  ports)  of  the 
UDC  circuit  in  a  as  a  function  of  frequency.  The  dotted  lines  represent  the 
couplings  between  different  ports.  The  solid  and  the  dashed  arrows 
represent  different  frequencies  and  respective  conjugates.  The  resonance 
line  shapes  of  the  two  spatially  distinct  ports  A  and  B  are  centred  at 
a>/\  =  mb  =  1  AAa.bCa.b-  ^ere  we  show  the  case  when  the  incoming  signal 
at  a>o  is  resonant  with  the  centre  frequency  (too  =oia.b)-  The  two 
sidebands  generated  by  the  UDC  on  channel  C  are  detuned  from  the  pump 
at  u)c  by  equal  amounts. 


The  main  building  block  of  our  design  is  a  reversible 
in-phase/quadrature  modulator  capable  of  carrying  out  noise¬ 
less  frequency  up-  and  down-conversion.  A  convenient  analyt¬ 
ical  model  capturing  the  fundamental  properties  of  the  device 
is  shown  in  Fig.  2a.  The  device  comprises  two  low-frequency 
inductor-capacitor  (LC)  resonators  (addressed  by  two  semi-infinite 
transmission  lines  A  and  B)  coupled  to  a  high-frequency  resonator 
(addressed  by  the  transmission  line  C)  through  time-varying  cou¬ 
plings  Mi, M2  that  emulate  the  role  of  the  pump  drive  in  active 
nonlinear  devices  and  transfer  energy  from  the  tone  at  a)c  to  the 
signal  modes  propagating  on  the  transmission  lines.  It  operates 
in  a  manner  analogous  to  the  in-phase/quadrature  modulation 
schemes  routinely  used  in  radiofrequency  communication  systems 
and  microwave  pulse  engineering  (hence  the  name)  and  converts 
two  orthogonal  spatial  modes  travelling  on  two  distinct  spatial  ports 
(A,  B)  at  the  same  frequency  (here  a> 0 )  into  two  orthogonal  temporal 
modes  travelling  on  the  same  spatial  line  (C)  at  different  frequencies 
(&>+ ,  ft )_).  However,  in  view  of  the  reversible  and  noiseless  frequency 
conversion  carried  out  by  this  device  (Fig.  2a),  we  will  henceforth 
refer  to  it  as  the  up/down-converter  (UDC).  In  practice,  such  a  de¬ 
vice  can  be  implemented  on-chip  using  two  ring  modulators  based 
on  Josephson  junctions,  along  the  lines  of  the  recently  demonstrated 
experiment  with  a  Josephson  parametric  converter17'18 . 

The  complete  design  for  the  active  circulator  (Fig.  3a)  consists  of 
a  UDC  functioning  as  a  frequency  up-converter,  a  phase  shifter  and 
a  second  UDC  functioning  as  a  frequency  down-converter. 

A  concise  representation  of  the  dynamics  at  each  of  the  three 
stages  in  the  cascade  is  provided  by  the  scattering  matrix  S 
that  relates  the  outgoing  wave  amplitudes  to  the  incoming  wave 
amplitudes  as  seen  from  the  ports  of  a  network.  We  begin  by 
deriving  the  scattering  matrix  of  the  UDC  stage  and  obtain  (see  the 
equations  (5)  and  (6)  of  Methods  section) 
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Here  a  and  b  denote  the  (reduced)  amplitudes  or  the  annihilation 
operators  for  the  waves  travelling  on  left  and  right  transmission 
lines  of  the  UDC  respectively.  This  notation  will  be  consistently 
used  to  denote  the  left-  (a)  and  right-  (b)  propagating  modes  at 
each  stage  of  the  device  hereafter  (see  Fig.  3a).  These  satisfy  bosonic 
commutation  relations  of  the  form19 


on  low-noise  microwave  parametric  amplifiers,  which,  with  the 
exception  of  designs  based  on  the  current-biased  d.c.  supercon¬ 
ducting  quantum  interference  device  (SQUID),  are  also  operated 
in  a  reflection  mode  with  both  the  input  and  output  signals 
collected  on  the  same  spatial  channel15'16.  However,  circulators 
(and  isolators)  routinely  use  bulk  components  made  from  ferrites 
to  achieve  non-reciprocal  phase  shifts  (Fig.  lc)  through  Faraday 
rotation,  making  them  unamenable  to  chip  fabrication.  Moreover, 
to  bias  the  magnetic  field  in  the  ferrite,  most  of  these  devices  use  a 
permanent  magnet,  which  may  channel  flux  into  the  superconduct¬ 
ing  device  under  test. 

Here,  we  present  the  full  analysis  of  a  model  for  a  four-port 
circulator  based  on  parametric  active  devices  with  no  magnetic 
components.  In  active  devices  the  energy  source — provided  by  the 
pump — acts  as  the  external  ‘bias’  field  and  sets  the  reference  phase 
for  the  system,  in  analogy  with  the  role  played  by  the  magnetic 
field  in  a  Faraday  medium.  We  exploit  this  effect  in  a  cascade  of 
active  devices  with  pump  phases  at  each  stage  shifted  appropriately 
to  obtain  non-reciprocal  transmission. 


[ a “/out,  a/'in/out]  =  S(cOi  -  a>j)  (2) 

In  writing  equations  (1)  and  (2),  we  have  set  a0  =  a[co0],a+  = 
a[ft>+],a_  =  «[&>_]  (see  Fig.  2b).  Similarly,  the  reflection  coefficients 
at  various  ports  are  denoted  by  r0,  r+  and  r_.  The  cross  reflection 
between  the  low-frequency  signal  ports  is  denoted  by  q0.  The 
transmission  coefficients  are  written  as  t  (transmission  without 
conjugation)  and  s  (transmission  with  conjugation)  with  subscripts 
(u,  d)  indicating  the  up-conversion  and  down-conversion,  respec¬ 
tively.  These  scattering  coefficients  are  a  function  of  frequency  in 
general  (see  Supplementary  Information  for  a  detailed  derivation). 
It  is  useful  to  note  that  the  phase  of  the  carrier,  denoted  by  (j> ,  affects 
only  the  transmitted  amplitudes  and  rotates  the  two  sidebands  in 
opposite  directions  as  can  be  seen  from  the  corresponding  scattering 
coefficients  5  and  t  in  equation  (1).  The  invariance  of  reflection 
amplitudes  to  the  phase  of  the  coupling  will  be  important  in 
understanding  total  reflections  of  the  cascade,  as  we  describe  later. 

Furthermore,  we  note  that  the  matrix  obtained  in  equation  (1) 
is  non-unitary,  that  is  SfS  /  1,  which  implies  non-conservation 
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Figure  3  |  Description  of  the  active  circulator,  a,  Circuit  schematic  of  the  active  circulator  design.  The  first  UDC  stage  acts  as  a  frequency  up-converter 
(UC;  also  indicated  by  a  gradation  in  the  colour  of  the  relevant  square)  with  a  parametric  coupling  modulated  at  the  carrier  frequency 
a>c  =  a>+  —  aio  =  ai-  +010  and  a  phase  4>l  =  0.  This  is  followed  by  a  phase  shifter  (PS)  that  phase  shifts  both  the  sidebands  by  n/2,  in  opposite  directions. 
They  are  then  demodulated  by  the  final  UDC  stage  acting  as  a  frequency  down-converter  (DC),  with  the  pump  phase  <p r  =  it/4,  b,  Forward  (green)  and 
backward  (red)  propagation  diagrams  calculated  using  the  transfer  matrix  method  for  a  with  an  appropriate  choice  of  detuning  (<5-i-  =  1  /-s/2)  and  coupling 
strengths  (q>l  =  “r  =  2-3/4)  for  maximum  isolation,  c,  Representation  of  the  device  operation  using  modulation  ellipses  at  each  stage  in  the  cascade.  The 
top  panel  shows  the  forward  propagation  for  the  case  when  two  distinct  signals  enter  ports  1  and  2  respectively  (do  =1,oo  =  0.5e1'';  tj  =  n/24)  and  the 
bottom  panel  shows  the  backward  propagation  dynamics  when  the  same  signals  enter  ports  3  and  4  respectively  (bo  =  1,bo  =  0.5e"U.  The  relative  phase 
and  amplitudes  are  chosen  to  represent  the  most  general  case  of  two  input  signals  that  differ  in  both  amplitude  and  phase.  The  relative  phase  difference 
between  the  two  signals  is  encoded  as  the  tilt  of  the  modulation  ellipse  in  the  in-phase/quadrature  plane  and  their  average  phase  is  represented  as  the 
colour  along  the  perimeter  of  the  ellipse,  with  yellow  indicating  zero  phase  (also  see  Supplementary  Information  for  more  details). 


of  photon  number  as  is  natural  for  an  active  device.  The  matrix 
recovers  its  unitary  form  as  we  turn  off  the  couplings  Mi,M2 
responsible  for  energy  transfer  between  the  pump  and  the  signal 
modes.  The  full  8x8  matrix  (Supplementary  Equation  SI 7),  de¬ 
scribing  the  device  operation  for  all  modes  and  their  respective  con¬ 
jugates,  fulfils  the  fundamental  requirement  of  symplecticity17,20. 

We  can  similarly  describe  the  action  of  the  frequency- 
independent  phase-shifting  stage  using  a  scattering  matrix  of 
the  form 


/  a+out  \ 
n_tout 
b+out 

V  b_fout  / 


/  0  0 

0  0 

e-,e  0 
V  0  eie 


~ie  0  \ 

0  e'e 

0  0 

0  0  / 


/«+m\ 

aJm 

b+in 

\bJmJ 


(3) 


The  scattering  matrix  of  the  whole  device  is  then  obtained  using  the 
transfer  matrix  formulation,  as  described  in  the  Methods  section 
(also  see  Supplementary  Information). 

In  our  analysis  we  consider  the  operation  at  resonance,  that  is, 
when  the  input  signal  frequency  coincides  with  the  band  centre  of 
the  input  resonators.  Setting  the  phase  of  the  pump  at  the  first  UDC 
stage  4> l  =  0  for  calculational  simplicity,  we  observe  a  transmission 
resonance  for  8  =  ±Jt/2  (phase  rotation  by  the  phase  shifting  stage), 
(j>R  =  j[/4  (phase  of  the  pump  at  the  second  UDC  stage),  =  l/v/2 
(detuning  of  the  sidebands  from  the  pump  in  units  of  linewidth, 
that  is,  half-width  at  full-maximum  of  the  resonance  line  shape)  of 
the  high-frequency  resonator)  and  aR  =  aR  =  M0/^/La,bTc  =  2~3/4 


(strength  of  the  parametric  coupling) .  For  this  choice  of  parameters, 
we  obtain  the  scattering  matrix  of  the  complete  device  as 

/  0  0  0  i  \ 

c  =  0  0  -i  0  (A\ 

itotal  i  0  0  0  1  J 

VO  i  0  0/ 

This  is  the  matrix  of  a  perfect  four-port  circulator.  The  analogy 
between  a  conventional  circulator  and  the  active  circulator  design 
proposed  in  this  Letter  is  made  apparent  from  the  respective 
wave  propagation  diagrams  in  Figs  lc  and  3b  (see  Supplementary 
Information  for  details  on  the  calculation  of  coefficients  on  different 
arms  in  Fig.  3b).  Nonetheless  there  are  important  differences 
between  the  two  designs  despite  the  identity  of  the  final  S  matrix. 
The  coefficients  on  the  forward-  (green)  and  backward-  (red) 
propagating  arms  of  the  active  circulator  design  (Fig.  3b)  involve 
deamplification  followed  by  amplification,  unlike  the  passive 
splitters  (90°  or  180°  hybrids)  employed  in  Faraday  rotation 
schemes.  This  can  be  observed  by  squaring  the  amplitudes  on  each 
of  the  two  arms  originating  from  (or  terminating  into)  a  port  and 
calculating  the  net  power  output,  for  each  isolated  UDC  stage.  It  is 
straightforward  to  observe  that,  unlike  the  case  of  Fig.  lc,  they  do 
not  add  up  to  unity.  The  overall  transmission,  however,  is  still  unity 
owing  to  an  exact  cancellation  of  the  larger  and  smaller  than  unity 
gains.  The  wave  propagation  diagrams  in  Fig.  3b  reveal  another 
important  difference  of  this  design  from  that  of  a  conventional 
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Figure  4  |  Variation  of  the  difference  between  forward  and  backward  transmission  coefficients  (IS31I  —  |si3 1 ).  a-c,  Asymmetry  in  transmission  calculated 
as  a  function  of  strength  of  the  coupling  a  and  phase  rotation  9  carried  out  by  the  second  phase-shifting  stage  (a),  detuning  <5-t  of  the  sidebands  from  the 
pump  and  phase  rotation  9  (b),  and  detuning  <5±  and  the  phase  of  the  pump  at  second  UDC  stage  0 r  (c).  The  points  of  maxima  correspond  to  the  ideal 
values  reported  in  the  text.  The  plot  in  b  also  shows  the  periodicity  of  the  response  of  the  device  as  a  function  of  9.  In  c,  the  variation  with  respect  to  the 
pump  phase  shows  the  reversal  of  transmission  characteristics  with  0r  — >•  — 0r.  As  in  b,  the  response  is  periodic  in  0r  with  a  period  equal  to  it.  It  can  be 
seen  that  the  design  continues  to  work  for  moderate  deviations  from  the  preferred  phase  angle  9  =  it/2,  coupling  ai  =  ctr  =  2-3/4,  detuning  =  1/V2  and 
pump  phase  0r  =  jt/4  (values  indicated  with  dashed  arrows  along  the  axis). 


circulator.  The  non-reciprocal  action  of  the  active  circulator  is  not 
based  on  any  non-reciprocal  phase  shifters;  instead  it  relies  on  the 
active  stages  used  for  frequency  up-  and  down-conversion.  The 
phase-matching  condition  in  the  forward  direction  is  met  by  tuning 
the  phase  of  the  coupling  at  the  input  and  output  UDC  stages.  In  the 
reverse  direction,  the  phase  mismatch  leads  to  unity  transmission  in 
the  spatially  orthogonal  port  instead,  leading  to  complete  isolation 
between  the  incident  signal  port  and  its  corresponding  output  port. 
Figure  3c  shows  a  convenient  method  to  visualize  this  circulator 
action  geometrically  by  mapping  the  device  dynamics  at  different 
stages  using  a  modulation  ellipse  representation  (see  the  Methods 
section  and  Supplementary  Fig.  S2  for  details). 

Furthermore,  as  seen  from  equation  (1),  the  reflection  coeffi¬ 
cients  at  the  UDC  stages  are  non-zero  for  all  modes.  Flowever, 
for  the  whole  cascade,  the  total  reflection  is  identically  zero  at 
every  port  (su  =  0  for  all  i  in  equation  (4)).  This  remarkable 
cancellation  of  total  reflections  for  the  cascade  can  be  understood 
in  analogy  with  a  Fabry-Perot  resonance  where  a  cavity  flanked 
by  two  identical  reflecting  mirrors  exhibits  unity  transmission  at 
resonance.  The  total  phase  shift  between  the  active  ‘mirrors’  in  our 
device:  (jc/2)0+  —  (— jc/2)n  t  =  Jt,  is  akin  to  the  resonance  condition 
when  a  half-wavelength  of  the  incident  radiation  equals  the  length 
of  the  Fabry-Perot  cavity.  Also  the  reflections  at  the  two  UDC  stages 
are  identical  (as  the  reflection  coefficients  are  independent  of  the 
phase  angle  <0,  equation  (1)),  fulfilling  the  second  condition  for  the 
transmission  resonance  and  net  cancellation  of  reflections.  A  point 
of  distinction  between  the  two  pictures  is  that  the  symmetry  of  the 
active  circulator  design  is  described  by  a  subgroup  of  the  SU( 4) 
group  (the  group  formed  by  4  x  4  complex  matrices  of  unit  deter¬ 
minant,  equation  (4)),  and  not  by  the  SU( 2)  group  that  describes 
passive  lossless  two-port  devices  such  as  the  Fabry-Perot  resonator. 

We  note  a  recent  theoretical  paper21  that  also  showed  the 
existence  of  noiseless  non-reciprocity  in  a  nonlinear  circuit  in¬ 
volving  microwave  resonators  coupled  through  Josephson  junc¬ 
tions.  Although  this  proposal  is  based  on  a  passive  rather  than 
active  Josephson  circuit,  unlike  our  proposal,  it  involves  very  small 
Josephson  junctions  susceptible  to  offset  charges.  Also,  owing  to 
its  low  characteristic  energy,  it  handles  qubit  readout  signals  with 
a  lower  dynamic  range  than  our  proposed  device.  Furthermore, 
the  dependence  of  the  action  of  our  active  circulator  design  on 
different  parameters  (Fig.  4)  shows  that  the  isolation  achieved  is 
robust  to  reasonable  deviations  of  parameters  from  their  ideal 
values;  thus,  the  active  circulator  design  holds  promise  for  use  in 
practical  circuits.  Another  interesting  feature  of  our  device  is  the 


reversal  of  its  transmission  characteristics  with  the  phase  of  the 
pumps  (0> l.r  !->■  — <0l,r)  (Fig.  4c).  In  the  classic  circulators  based  on 
passive  Faraday  rotation,  this  can  be  accomplished  by  changing  the 
polarity  of  the  magnetic  bias  field.  Thus,  the  relative  shift  of  the 
pump  phase  in  an  active  device  indeed  plays  a  role  equivalent  to 
the  magnetic  field  in  a  Faraday  medium. 

Our  scheme  for  non-reciprocal  wave  propagation  relies  on 
the  directionality  imposed  by  shifting  the  phases  of  the  pump 
drives  in  a  cascade  of  parametrically  active  stages.  The  proposed 
design  operates  noiselessly  as  the  complete  device  has  unity  gain 
and  consists  of  purely  dispersive  components  with  no  dissipation, 
making  it  attractive  for  quantum  information  applications  using 
superconducting  circuits22.  As  the  proposed  protocol  is  generic  in 
nature,  it  should  be  easily  adaptable  to  optical  frequencies  where  it 
can  complement  the  recently  proposed  designs  of  non-reciprocal 
light  propagation  based  on  dynamical  modulation  of  the  refractive 
index  of  photonic  structures23  and  the  use  of  a  surface  waveguide 
on  photonic  crystals24.  In  addition  to  providing  a  design  for 
practical  on-chip  non-reciprocal  devices,  the  treatment  described 
in  this  Letter  may  also  give  theoretical  insights  into  the  inherently 
directional  dynamics  of  devices  such  as  the  microwave  d.c.  SQUID 
(ref.  25),  when  additional  active  stages  are  included  in  the  chain. 
This  would  be  useful  in  tackling  unanswered  questions  pertaining 
to  the  quantum  noise  of  d.c.  SQUID  amplifiers26. 

Methods 

The  scattering  matrix  for  the  UDC  is  calculated  from  the  network  impedance 
matrix  Z,  as  seen  from  its  ports,  using  the  identity 

S  =  (Z  +  Z„)-1x(Z-Z„)  (5) 

where 

Zo  =  diag(ZA,ZB,Zc,Zc)  (6) 

Here,  the  symbol  x  denotes  matrix  multiplication  and  Za  =  Zb  ,  Zc  are  the 
characteristic  impedances  of  the  semi-infinite  transmission  lines  serving  as  low- 
and  high-frequency  ports  respectively. 

For  each  stage  of  the  cascade,  we  now  transform  from  the  scattering  matrix, 
which  provides  a  port-neutral  relation  between  various  input  and  output  ports,  to 
a  port-specific  transfer  matrix: 


The  transfer  matrix  representation6  is  especially  useful  for  analysing  the 
proposed  design,  because  it  is  straightforward  to  calculate  the  total  transfer 
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matrix  of  the  device  by  multiplying  the  respective  transfer  matrices  of  different 
stages  in  the  cascade, 

Total  =  TDCr  X  Tps  x  7\jcl  (8) 

Here  the  subscripts  L,  R  index  the  left-hand  upconversion  (UC)  and  right-hand 
downconversion  (DC)  stage  and  the  subscript  PS  refers  to  the  intermediate  phase 
shifter  (Fig.  3a).  The  scattering  matrix  of  the  whole  device  is  obtained  from 
Total  using  the  inverse  of  the  transformation  in  equation  (7).  For  the  purpose 
of  evaluating  Total  from  equation  (8),  Tuc,  T> c  and  TPS  are  calculated  from 
equations  (1)  and  (3)  respectively,  while  noting  that 

T> c  =  .F-1  x  Tjc-1  x  F 

=  PxToc-1xF,  (F~'=F)  (9) 

In  equation  (9),  F  =  crx  ®  T  (ox  is  the  two-dimensional  Pauli  spin  matrix  and 
I2  is  the  two-dimensional  unity  matrix).  This  matrix  F  is  required  to  flip  the 
indices,  thus  maintaining  consistency  in  labelling  the  ‘in’  and  ‘out’  amplitudes 
along  a  given  direction  of  propagation  (see  Supplementary  Information  for 
additional  details). 

The  modulation  ellipses  in  Fig.  3c  are  calculated  using  the  transfer  matrices, 
with  coefficients  indicated  in  Fig.  3b.  This  representation  is  inspired  by  the 
polarization  ellipse  used  to  represent  the  state  of  polarization  of  an  electromagnetic 
wave  (linear,  circular  or  elliptical),  which  involves  recording  the  trajectory  traced 
out  by  the  tip  of  the  polarization  vector  of  light  (defined  by  the  instantaneous 
direction  of  the  electric  field  vector  E)  in  a  plane  perpendicular  to  the  direction  of 
propagation.  Equivalently,  a  two-dimensional  representation  of  the  components  Ex 
and  EY  of  the  electric  field  in  the  complex  plane  can  be  used  to  obtain  a  geometric 
description  of  the  polarization  of  the  light  wave.  In  the  case  of  a  modulation  ellipse 
representation  of  Fig.  3c,  we  extend  this  idea  to  map  two  distinct  orthogonal 
modes  (x,y)  at  each  stage  of  the  proposed  device  (spatial:  (x  =  ao ,y  =  af  )  or 
temporal:  (x  =  a+,y  =  a_))  as  an  ellipse  in  the  plane  defined  by  the  coordinates 
I  =  Re[x  +y],  Q  =  Im[x  —  y*].  This  exercise  shows  that  the  final  ellipses  obtained 
at  the  output  in  the  cases  of  forward  and  backward  propagation  through  the  device 
are  rotated  by  90°  with  respect  to  each  other,  indicating  that  in  the  case  of  reverse 
propagation  the  orthogonal  spatial  port,  relative  to  the  forward  propagation, 
receives  the  transmitted  energy,  leading  to  a  circulator  action. 
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